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Abstract: It was recently pointed out that in super symmetric scenarios with gravitino 
dark matter and bilinear R-parity violation, gravitinos with masses below My/ typically 
decay with a sizable branching ratio into the 3-body final states W*£ and Z*v. In this 
paper we study the indirect detection signatures of gravitino dark matter including such 
final states. First, we obtain the gamma ray spectrum from gravitino decays, which features 
a monochromatic contribution from the decay into 71/ and a continuum contribution from 
the three-body decays. After studying its dependence on supersymmetric parameters, we 
compute the expected gamma ray fluxes and derive new constraints, from recent FERMI 
data, on the R-parity breaking parameter and on the gravitino lifetime. Indirect detection 
via antimatter searches, a new possibility brought about by the three-body final states, 
is also analyzed. For models compatible with the gamma ray observations, the positron 
signal is found to be negligible whereas the antiproton one can be significant. 
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1. Introduction 

Supersymmetric scenarios are considered to be the most promising extensions of the Stan- 
dard Model. They are motivated by the hierarchy problem, by the unification of the gauge 
couplings, and by the existence of dark matter. Recent observations indicate that dark 
matter accounts for about 25% of the energy-density of the Universe Q, a fact that cannot 
be explained within the standard model of particle physics. Supersymmetric models, on 
the other hand, contain at least two viable dark matter candidates: the neutralino and 
the gravitino. Both give rise to an interesting phenomenology and have been extensively 
studied in the literature. In this paper we focus on models with gravitino dark matter. 

In supersymmetric models, R-parity conservation was long believed to be a prerequisite 
for supersymmetric dark matter, for it guarantees the stability of the lightest supersymmet- 
ric particle (LSP) -the would-be dark matter candidate. It was pointed out in ||], however, 
that if the gravitino -the superpartner of the graviton that arises in local supersymmetric 
theories- is the LSP, it can be a suitable dark matter candidate in such scenarios. In- 
deed, even though the gravitino is unstable due to R-parity breaking, its lifetime may be 
much longer than the age of the Universe. Such a long lifetime is the result of the grav- 
itino weak interactions, which are suppressed by the Planck scale and by small R-parity 



violating couplings. A scenario with gravitino dark matter and R-parity violation is also 
favoured by thermal leptogenesis, as it alleviates the tension between the large reheating 
temperature required by leptogenesis and the constraints from Big-Bang Nucleosynthesis 
Q. In R-parity violating models, therefore, the primordial gravitino produced in the early 
Universe is a viable and well-motivated dark matter candidate. 

A salient feature of such models is that, unlike their R-parity conserving counterparts, 
gravitino dark matter can be indirectly detected. In fact, the small fraction of gravitinos 
that have decayed until today constitutes a source of high energy cosmic rays [||, |H| that 
could be observed in present and future experiments. Gravitinos may be indirectly detected 
via gamma rays [Q, y], neutrinos [|7J or antimatter searches ||. Present data from FERMI, 
for instance, already constrains the parameter space of these models j|]. 

The indirect detection signatures of gravitino dark matter strongly depend on the 
gravitino lifetime and on its branching ratios. It was recently demonstrated that, in bilinear 
R-parity breaking models, gravitinos with masses below Myy not only decay into 7^, as 
previously believed, but also into the three-body final states W*£ (— > Ylfff'ty an d Z* v 
(— > Ylf i f v ) 01 • Because the branching ratio into these three-body final states is typically 
sizable, they considerably alter the gravitino lifetime and the indirect detection signatures 
of gravitino dark matter. In this paper we revisit the prospects for the indirect detection 
of gravitino dark matter in view of these new decay modes of the gravitino. 

In the next section we introduce the model and explain the main assumptions of our 
work. Then, we briefly revisited the dominant decay modes of gravitino dark matter. In 
section |4|, the expected gamma ray flux from gravitino decay is studied, in particular its 
dependence on gaugino and gravitino masses. Then, we use that flux and recent Fermi data 
to derive constraints on the R-parity breaking parameter and on the gravitino lifetime. 
In section ^ the indirect detection of gravitino dark matter via antimatter searches is 
considered. After deriving the spectra of positrons and antiprotons, we obtain the expected 
fluxes at earth and compare them with present data. For models compatible with the 
gamma ray observations, the positron signal is predicted to be negligible whereas the 
antiproton one can be significant. 

2. Bilinear R-parity violating models 

We work in the framework of supersymmetric models with bilinear R-parity violation and 
gravitino dark matter, such as those considered in y, || g]. That is, we assume that the 
gravitino is the lightest supersymmetric particle and that it accounts for the observed dark 
matter density of the Universe, Q^h 2 ~ 0.11. In these models, the gravitino is unstable, 
due to R-parity breaking, and may decay into standard model particles via lepton number 
violating interactions. 

In bilinear R-parity violating models, the superpotential is given by 

W = Wmssm + ^L t H u , (2.1) 

where Wmssm is the usual (R-parity conserving) superpotential of the MSSM, Li (i = 
e, //, t) is the lepton doublet, H u is the up-type Higgs doublet, and /ij are R-parity violating 



couplings. To study these models, it is convenient to redefine the Higgs and lepton doublets 
so as to eliminate the bilinear R-parity violating interactions from the superpotential. In 
that basis, which is the one we consider in the following, the mixing terms between the 
Higgsino and the lepton doublet are eliminated from the fermion mass matrix, and the 
soft-breaking Lagrangian contains the following terms 

£soft = 4^ t SM + BiLiH u + m\ H LiH* d + h.c. (2.2) 

Here Hj is the down type Higgs boson doublet, and Bi,m 2 f are R-parity violating 



parameters. Even though the original superpotential, equation (2.1), did not contain any 
trilinear R-parity violating terms, the redefinition of the fields that eliminates the bilinear 
ones induces A and A' terms in the superpotential. They are given by 

Wrpv = XijkLkLtEj + X ijk L k QiDj , (2.3) 

where Lj and Qi are the lepton and quark doublets, whereas Ef and Df are the lepton and 
down-quark singlets. The couplings A and A' are related to the original couplings of the 
superpotential by 

\ ijk = ^\f J and A^ = ^A£, (2.4) 

fJ, (J, 

where \i is the usual MSSM parameter and A e and A are respectively the charged lepton 
and down-quark Yukawa couplings. As we will see, these trilinear terms play essentially 
no role in our discussion. 

The two R-parity violating terms present in £ so ft give rise to non-zero vacuum expec- 
tation values for the sneutrino fields. They can be written as 

Bi sin (3 + m 2 ~ cos (3 

(h) = r^ v ' ( 2 - 5 ) 

where v = 174 GeV, tan/3 = (H®)/(H®), and m^ i is the sneutrino mass. It is useful to 
define the dimensionless parameters £$, % = e,fi, r, as 

6 = ^- (2-6) 

v 

Because the R-parity violating couplings are expected to be largest for the third generation, 
we will assume, following [7j, that the sneutrino acquires a vev only along the v T direction. 
In other words, we suppose that £ T ^> £^,£ e - I n this setup, all R-parity violating effects, 
including the decay of the gravitino and non-zero neutrino masses, are controlled by £ r . 

Since neutrino masses may arise from other sources, like the seesaw mechanism, we 
will only require that the contribution to neutrino masses from R-parity violation should 
not exceed the value obtained from present data. This condition implies that £ r < 10 . 
In addition, a lower bound on £ T is obtained from successful big bang nucleosynthesis: 
£ r > 10~ n ||. A priori, therefore, £ T can vary over four orders of magnitude. In section [| 
we will use the gamma ray flux from gravitino decays and recent Fermi data to derive new 
constraints on £ r . 



3. Gravitino decays 

In this section we briefly revisit the decays of the gravitino in bilinear R-parity violating 
scenarios. We assume throughout this paper that m^ < Myy, for it is only in such mass 
range that the three-body final states modify the gravitino lifetime and branching ratios. 



The gravitino decay rates can be calculated directly from the interaction Lagrangian [11]. 
It turns out that gravitinos with masses below My/ can decay into the two-body final state 
7^ T or into the three-body final states W*t and Z*u T . The decays induced by the trilinear 



couplings, equation (2.4), are suppressed by the Yukawa couplings and by sfermion masses 
and turn out to be negligible compared to these (see e.g. |J), so we will not considered 
them in the following. For the two-body final states 71V, the gravitino decay width is given 

by § § 



T(G -> -yv r 



& m G irr 12 



(3.1) 



where U~^ is the photino-zino mixing parameter 



written as ]|, 12] 



\U~y\ 



M Z (M 2 



. To a good approximation it can be 
Mi)s w c w 



(M lC 2 + M 2 sl)(M lS 2 w + M 2 c 2 w ) 



(3.2) 



where sy/ and cy/ denote the sine and cosine of the weak mixing angle, and Mi,M 2 are 
the U(l) and SU(2) gaugino masses. For simplicity, we will assume in this paper that 
gaugino masses are universal at the GUT scale, so that M 2 ~ 1.9Mi at the electroweak 
scale, and we will consider Mi as a free parameter of the model. In any case, our results 
do not depend strongly on this universality assumption. Notice, from the above equations, 
that the mixing parameters and the decay rate into jv T decrease with gaugino masses. 




a 




Figure 1: The two diagrams that contribute to gravitino decay into W*t 



The gravitino decay rates into the three-body final states W*t and Z*v T were first 
computed in |l0|. Each of these processes receives contributions from two different dia- 
grams, as illustrated in figure [l] for the decay into W*t. In the appendix we provide the 
analytical expressions for both decay rates. A crucial feature of these decay rates is that 
they contain terms independent of the wino-wino and zino-zino mixing parameters that 



dominate the total rate for large gaugino masses. The decay into W*t is always more im- 
portant than that into Z*u T and they typically give a significant correction to the gravitino 
decay width, particularly for large gaugino masses. Figure ^ shows the gravitino branching 
ratio into the three-body final states, BR(G — > W / *r)+BR(G — >■ Z*v T ), as a function of the 
gravitino mass for different values of M\. As clearly seen in the figure, the three-body final 
states give an important contribution to the decays of the gravitino. They even dominate 
gravitino decays over a wide range of gaugino and gravitino masses. 
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Figure 2: The gravitino branching ratio into three-body final states (BR(G — > W*t)+BR(G — > 
Z*v T )) as a function of the gravitino mass for different values of M\. These branching ratios are 
independent of £ T . 

The total decay width of a gravitino with mg < Myy will, then, be given by 

r tot (G) = 2(T(G ->■ 7 i/ T ) + r(G -)- r+IF"*) + T(G -► u T Z*)), (3.3) 



where the factor two takes into account the charge conjugated final states. This decay 
width is proportional to the R-parity breaking parameter £j!. The gravitino lifetime, tq, is 
simply the inverse of the decay width, l/T to t(G). 

Because the indirect detection signatures of gravitino dark matter depend on the grav- 
itino lifetime and its branching ratios, and they are both strongly affected by the three-body 
final states W*t and Z*v T , it is necessary to revisit the prospects for indirect detection of 
gravitinos in view of these new decay modes. At a qualitative level, it is easy to understand 
the novel features induced by the decays into W*t and Z*v T . Regarding gamma rays, the 
implications are twofold. On the one hand, one expects a suppression in the gamma ray line 
from the two-body decay, due to the smaller branching into 71/,- associated to a given grav- 
itino lifetime. On the other hand, there will also be a new continuum contribution entirely 
due to the three-body final states. Regarding antimatter signals, the three-body gravitino 
decays give rise to a non-zero positron and antiproton flux for mg < M\y, opening up the 
possibility of constraining or discovering gravitino dark matter via antimatter searches. 



In the next sections we will study quantitatively the implications of the three-body final 
states for the indirect detection of gravitino dark matter. 

4. Gamma rays from gravitino decay 

Since the gravitino typically decays with sizable branching ratios into both ^u T and W*t + 
Z*v T , the gamma ray spectrum from gravitino decays features a monochromatic contri- 
bution from the two-body final state and a continuum contribution from the three-body 
final states. The relative weight between them is determined by the decay branching ra- 
tios, which depend on mg and M\, whereas the overall normalization of the spectrum is 
controlled by the R-parity breaking parameter £ T . In this section we study in detail the 
gamma ray spectrum from gravitino decay and we determine the constraints that can be 
imposed on scenarios with gravitino dark matter from recent Fermi data. 



4.1 The gamma ray spectrum 

The contribution of the two-body final state to the photon spectrum, which is simply 
a delta function at E = mg/2, has already been considered in previous works -see e.g. 
|I|, [|. The contribution of the three-body final states W*t and Z*v T , on the other hand, 
is computed here for the first time. The photons in this case are produced after the decay 
and fragmentation of the final states and come mainly from r bremsstrahlung and neutral 
pion decays. To obtain the spectra we have used the event generator PYTHIA [13|. The 
gamma ray spectrum from G decay, dN/dE, features a monochromatic and a continuum 
contribution and can be expressed as 



dN ~ dN~< v ~ dN w * T ~ dN z * u 

BR(G -> 1V )^— + BR(G -> W*r)% + BR{G -> Z*i/ T )- 



dE. 



dE~, 



dE~, 



dE~, 



monochromatic 



continuum 



(4.1) 



Here, it is not necessary to include the inverse compton contribution to the gamma ray 



flux, as it is known to be negligible [14] for the range of gravitino masses we are considering. 
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Figure 3: The continuum gamma ray spec- 
trum from gravitino decay for mg = 70 GeV 
and different values of Mi . 
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Figure 4: The continuum gamma ray spec- 
trum from gravitino decay for Mi = 1 TeV 
and different values of the gravitino mass. 



To begin with, let us study the continuum contribution to the gamma ray spectrum 
from gravitino decay. Its dependence on gaugino and gravitino masses is illustrated in 
figures H and [|. Figure [| shows this contribution as a function of the energy for m^ = 
70GeV and different values of M\. Because of the higher branching ratio into W*t and 
Z*v T , the continuum contribution increases with gaugino masses. This increase is quite 
important for low gaugino masses (notice the difference between the M\ = 100 GeV and 
the M\ = 300 GeV lines) but becomes irrelevant above M\ ~ 1 TeV. The reason is that 
the branching into the three-body final states is already very close to 1 in that region, so 
the spectrum hardly changes for even higher values of M\. 

Figure |] illustrates the variation of the continuum contribution with the gravitino 
mass. In this figure M\ was set to 1 TeV, which results, for the gravitino masses we 
consider, in a three-body branching larger than 90% -see figure ^. We observe that the 
continuum contribution is sensitive to the gravitino mass at high energies and that it is 
almost independent of it for E < 1 GeV. 

Now that we have obtained the gamma ray spectra from gravitino decays, the next step 
is the computation of the expected flux at earth. Two different populations of gravitinos 
contribute to such flux: gravitinos from the Milky Way and gravitinos at cosmological 
distances. It turns out that the dominant contribution comes from gravitinos in the Galactic 
halo y, 14] so we will focus on those. The gamma ray flux at earth from gravitinos decaying 



in the Milky Way halo can be computed as 

# 7 (V>) 2 dN 1 



dE mg dEj 87rrg ji s 



p(l)dl, (4.2) 



where p(l) is the dark matter distribution in the halo, and the integral extends over the line 
of sight. Notice from this equation that the gamma ray flux from gravitino decay depends, 
on the particle physics side, only on the spectra, the gravitino mass and the gravitino 
lifetime, which is in turn determined by mg, £ T and M\. For p we use the Navarro, Frenk 
and White (NFW) profile |J], 

'M = / ivi i v > ( 4 - 3 ) 

r/r c (l +r/r c ) z 

with r c = 20 kpc and p c = 0.34 GeV cm -3 , which correspond to a local dark matter 
density po = p{r = 8.5 kpc) of 0.39 GeV cm~ 3 [16|. To facilitate the comparison with 



current observations, in our numerical analysis we will be concerned only with the high 
latitude region, \b\ > 10°. In contrast to the Galactic Center, this region has the benefit of 
not being that sensitive to the assumed halo profile. For that region we find that 



~d~E 



2 dN 1 ,, 



, 6 | >10 o m d dE 1 %KT d 



x 2.1 x 10-" GeV/Wstr. (4.4) 



Figures |5| and |6| show the dependence of the gamma ray flux from gravitino decay on mg 
and M\. The monochromatic contribution has been convoluted with a 10% gaussian energy 
resolution, comparable to that expected from FERMI. From these two figures we can see 
that the monochromatic and continuum contributions are clearly discernible. The former 
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dominates the flux at high energies, E ~ tti^/2, whereas the latter does it at somewhat 
lower energies. Figure || shows that the height of the monochromatic contribution is quite 
sensitive to the value of M\. The continuum contribution, on the other hand, is only slightly 
affected by M\. As expected from the gravitino branching ratios, the larger the gaugino 
masses the smaller the line signal and the larger the continuum. From figure ^ we see that 
m,Q not only determines the position of the line signal but also the height of the continuum 
contribution, being larger for heavier gravitinos. If such a flux were observed, m?; could be 
extracted from the position of the peak, while the gaugino mass Mi could be determined 
from the relative weight between the monochromatic and continuum contributions. 

In previous studies on the indirect detection of gravitino dark matter only the two-body 
decays of the gravitino were considered. That is, the branching ratio into ^v T was assumed 
to be equal to 1 for mg < My/ -see e.g. @, [|. It is useful, therefore, to compare the correct 
gamma ray flux expected from gravitino decays (two- and three-body decays included) with 
that obtained in previous works (only the two-body decay). Figure |7| shows the predicted 
gamma ray flux in both cases. For that figure we set mg = 70 GeV, M\ = ITeV, and 
tp; = 8 x 10 24 s. The correct flux features a much smaller monochromatic contribution and 
the entirely new continuum contribution. Clearly, to reliably assess the indirect detection 
prospects of gravitino dark matter via gamma rays, the three-body decays of the gravitino 
must necessarily be taken into account. 

Next, the gamma ray fluxes we just obtained are used in conjunction with recent 
FERMI data to derive new constraints on scenarios with gravitino dark matter and bilinear 
R-parity violation. 



4.2 Constraints on gravitino dark matter from FERMI data 

Recently, the Fermi collaboration published two results that can be used to constrain the 
parameter space of gravitino dark matter. On the one hand, they searched for photon lines 
with energies between 30 GeV and 200 GeV and derived a constraint on the annihilation or 
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Figure 7: A comparison between the correct gamma ray flux expected from a decaying gravitino 
(dashed line) and the corresponding flux obtained including 2-body final states only (dotted line). 
For this figure we set m^ — 70GeV, M\ = ITeVGeV, and Tg = 8x 10 24 s. The gamma ray 
background is also shown (solid line). 



decay rate of dark matter particles into monochromatic photons JIT]] 1 . Since the two-body 
decay of the gravitino produces a gamma ray line at mg/2, we can apply such constraints 
to our scenario. In addition, Fermi also measured the diffuse gamma ray emission for the 
Galactic latitude range \b\ > 10° and energies between 100 MeV and 100 GeV, and found 
it to be compatible with the theoretical expectation of the diffuse background obtained in 
models of cosmic ray propagation in our Galaxy |lq| . These measurements have already 
been used to constrain different models of decaying or annihilating dark matter, see e.g. 
J19| |. As we have seen, gravitino decays give an additional contribution to the gamma ray 
flux in that energy range, so we must make sure that the gravitino contribution is small 
enough to remain compatible with data. Here, we will use those two results to derive 
constraints on the value of £ T and on the gravitino lifetime for different values of Mi and 



m G- 



In figure |8] we show the exclusion regions, in the plane (rog,( T ), derived from the 
Fermi bounds on a gamma ray line from dark matter decay. Notice that, in our scenario, 
those bounds apply only to the range 60 GeV < mg < 80 GeV. To obtain the exclusion 



regions we simply rescaled the bound from Table I in 17] to take into account that only one 
photon is emitted in a gravitino decay. As the figure illustrates, the line constraint does 
not depend much on the gravitino mass -at least in the range of interest to our discussion. 
It does depend on gaugino masses. It is one order of magnitude weaker for M\ = 1 TeV 
than for M x = 100 GeV. 

We can also write these results as constraints on the gravitino lifetime rather than 
on £ r . Figure || shows exclusion regions in the plane (mg,rg) for different values of M\. 



1 A new analysis using two years of Fermi data has just been published 
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Figure 8: The constraint on £ T derived from the gamma ray line searches by Fermi. The region 
above the lines is already ruled out by Fermi data for the respective value of M\ . Notice that the 
FERMI gamma ray line searches |L7| do not constrain gravitino masses below 60 GeV. 
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Figure 9: The constraint on the gravitino lifetime derived from the gamma ray line searches by 
Fermi. The region below the lines is already ruled out by Fermi data for the respective value of 
Mi . Notice that the FERMI gamma ray line searches |17j do not constrain gravitino masses below 
60 GeV. 



The area below the lines is excluded by Fermi data for the respective value of M\. We see 
that the constraint on the gravitino lifetime lies approximately between rg > 10 29 s for 
M\ = 100 and r^ > 10 27 s for M\ = 1 TeV, with a slight dependence on the gravitino 
mass. As before, this bound on tq from Fermi line searches applies only to gravitino masses 
between 60 and 80 GeV. 

A bound valid over a wider mass range can be obtained from the isotropic diffuse Fermi 



10 



data. The spectrum of the extragalactic background measured by Fermi is compatible with 
a power law with index 7 = 2.41 ± 0.05 [18|. To obtain the exclusion constraints from the 



isotropic diffuse gamma-ray data [18[, we compare the expected signal plus background 
with Fermi data for each energy bin. A model is considered to be excluded if the signal 
plus background deviates from the data by more than 3<r in any energy bin. Exclusion 
zones in the plane (m^,^ T ) are obtained by increasing £ r until that condition is satisfied. 
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Figure 10: The maximum allowed values of £ T as a function of the gravitino mass for different 
values of M\. The region above the lines is already ruled out by Fermi data. 



Figure 10 shows the constraints on £ r derived from the diffuse Fermi data. It displays 
exclusion lines in the plane (mg, £ T ) for different values of M\. The region above the lines is 
excluded for the respective value of M\. The step- like features observed in the figure are the 
result of the line feature in the spectrum moving through the different FERMI energy bins. 
Notice that the constraint on £ r is more stringent at smaller gravitino masses and that it 
strongly depends on gaugino masses. It is, in fact, about one order of magnitude stronger 
for M\ = 100 GeV than for M\ = 1 TeV. This dependence on Mi clearly indicates that 
the monochromatic contribution is dominantly responsible for the constraint, in agreement 
with the fluxes we derived previously -see figures . Comparing figures 1C and || we notice 
that, for gravitino masses between 60 and 80 GeV, the strongest constraint on £ T comes 
from the line searches rather than the diffuse emission data. At smaller gravitino masses, 
where there are no bounds from line searches, the allowed values of £ T are determined by 
the diffuse Fermi data. 

The bound on £ T we have just obtained can be translated, via equation ( |3.3| ), into a 
constraint on the gravitino lifetime, as illustrated by Figure O. In it, the region below the 
lines is excluded by the diffuse Fermi data. Let us first emphasize that the dependence of 
these constraints on gaugino masses is entirely due to the three-body final states W*t and 
Z*v T . If only two-body final states were considered there would have been a unique line 
in this figure. From the figure we see that for M\ = 100 GeV the gravitino lifetime should 
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Figure 11: The minimum allowed values of tq as a function of the gravitino mass for different 
values of Mi . The region below the lines is already ruled out by Fermi data. 

be larger than about 2 x 10 28 s, with only a mild dependence on the gravitino mass. For 
Mi = 1 TeV, the constraint can be up to two orders of magnitude weaker, depending on 
the gravitino mass. As expected, at low gravitino masses the effect of the three-body final 
states is smaller (see figure 0), giving rise to a weaker dependence of the constraint with 
gaugino masses. 

To summarize, we have seen that the three-body decays of the gravitino considerably 
alter the gamma ray spectrum expected from G decays and the constraints that can be 
derived from gamma ray data. 



5. Antimatter searches 

Antimatter searches are another promising way to indirectly detect dark matter. The 
positrons and antiprotons produced in gravitino decays originate entirely in the W*t and 
Z*u T final states and had not, in consequence, been considered before in the literature. 
In this section, we obtain the positron and antiproton spectra from gravitino decay and 
compute their expected fluxes at earth. 

5.1 Positrons 

The positrons from gravitino decay into W*t and Z*v T originate in the decay of the r and 
in the fragmentation of the virtual gauge bosons, mainly via n + decay. They give rise to a 
continuum spectrum of positrons, which we obtain using the event generator PYTHIA [pi . 



Figure 12 shows the positron spectrum (before propagation) resulting from gravitino decay 
for top; = 70 GeV and different values of M\. At high energies, E 2 dN e +/dE is observed 
to increase with energy, indicating that dN e +/dE is a slowly decreasing function in that 
region -a result that can be attributed to the positrons from tau decay. As the gauginos 
become heavier, the branching ratio into three-body final states gets larger and so does the 
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positron flux. In fact, it is observed to be about four times larger for M\ = 1 TeV than 
for Mi = 100 GeV over the whole energy range. For M\ ~ 1 TeV, the branching into 
W*t + Z*v T is almost 1 and the positron flux cannot longer increase with gaugino masses. 
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Figure 12: The positron spectrum (before propagation) from gravitino decay for different values 
of Mi . For this figure we set mg = 70 GeV. 



Unlike gamma rays, which travel directly to us, positrons are affected by propagation 
effects in the Galaxy. A careful discussion of these effects and of the theoretical uncertain- 
ties involved was presented in p0| . The important point for us is that propagation effects 
can be taken into account by solving, for the positron number density / = dN e +/dE, the 
following stationary diffusiondoss equation: 





-K(E)V'f-—(b(E)f) = Q e+ 



(5.1) 



where K(E) = K (E/GeV) 5 is the diffusion coefficient and b(E) = E 2 /(GeV ■ t e ) is the 
energy loss coefficient, with t e = 10 16 s. They describe, respectively, transport through 
the turbulent magnetic field and energy losses due to synchrotron radiation and inverse 
compton scattering on CMB photons and on galactic starlight. For positrons originating 
in gravitino decays, the source term, Q e +, is given by 



Qe+(r,E) 



Top; *— ' 

Cr i 



dE 



(5.2) 



where the sum is over the different decay modes that can generate positrons. That is, over 
W*t and Z*v T . 



Equation (5T) is solved within a diffusive region with the shape of a flat cylinder. Its 
radius R is 20 kpc and its height is 2L in the z direction. On the surface of this cylinder, the 
positron density / vanishes. The solar system is located at (r, z) = (8.5 kpc, 0). Thanks to 
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these simplifying assumptions, the positron flux at earth can be written in a semi-analytical 



form (see e.g. |8|, 21]. In this framework, / will depend on the propagation parameters 
5, Kq and L, which are constrained by cosmic ray data. Here, we will simply adopt the 
so-called MED propagation model @ (5 = 0.70, K = 0.0112 kpc 2 /Myr,L = 4 kpc), 
which provides the best fit to the Boron-to-Carbon (B/C) ratio. 

To compare with experimental data, it is better to use, rather than the positron flux, 
the positron fraction: the ratio between the positron flux and the electron+positron flux. 
To compute the positron fraction, we need to know not only the positron flux from gravitino 
decays but also the background electron and positron fluxes. The interstellar background 
fluxes of electrons and positrons can be parameterized as |22] 

^- g (£) = ( -, 82 Z^L ) GeV-V-VV" 1 , (5.3) 



l + 0.224e 2 - 93 . 

*SW = ( l^OO^ + 24 - ^ 341 ) GeV-W VV- . (5.4) 

with e = Ej GeV. These expressions reproduce to an excellent approximation the results 
obtained for a GALPROP conventional model - Model in pq] — that is compatible with 
the gamma-ray diffuse emission spectrum measured by Fermi at intermediate Galactic 
latitudes. Finally, it must be taken into account that for energies smaller than about 
10 GeV the electron and positron fluxes at the top of the atmosphere (the measured ones) 
can differ, because of solar modulation effects, from the interstellar fluxes. In the force field 



approximation [ 24 ] , the fluxes at the top of the atmosphere are related to the intersterllar 



fluxes by the following equation: 

C ± OA (^toa) = ^^ S ± (^is) (5.5) 

where Ei$, -Etoa are the electron/positron fluxes at the heliospheric boundary and at 
the top of the Earth's atmosphere respectively. They are related with each other via 
the solar modulation parameter, (ftp, as Eis = -Etoa + <Pf- In our calculations we use 
4>F = 500 MV, which is characteristic of the minimum of solar cyclic activity. We now 
have all the ingredients required to compute the positron fraction in our scenario and to 
compare it with present data. 



Figure 13 shows the predicted positron fraction as a function of the energy for m^ = 
80 GeV, M\ = 1 TeV, and two different values of £ T . In addition, the PAMELA measure- 
ments ]2J| and the predicted background, from equations (|5.3| ) and G5,4| ), are also shown. 
As is well known, the PAMELA data indicates an excess of positrons over the expected 
background [25]. The inclusion of the three-body final states W*t and Z*v T does not help 



in resolving that discrepancy, as they affect the positron flux only at relatively low ener- 
gies, E < 30 GeV. From the figure we see that for £ T = 5 x 10 -8 , gravitino decays would 
produce a significant deviation in the positron fraction. Unfortunately, such large values of 
£ r are not compatible with gamma ray constraints -see figure [l(]. The largest value allowed 
by such data is about one order of magnitude smaller, £ T = 7.3 x 10 -9 . As illustrated in 



figure 13, such a small £ T does not give rise to an observable signal. A result that, as we 
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Figure 13: The positron fraction as a function of the energy for m^ = 80 GeV and M\ = ITeV. 
The interstellar background and the PAMELA data are also shown. Two different values of £ T are 
illustrated. For £ r = 5 x 10 -8 (dashed gray line) the positron fraction is clearly distinguishable from 
the background but such value of £ r is not compatible with gamma ray data. The solid line shows 
the positron fraction computed for the largest value of £ T allowed by the gamma ray constraints 
(7.3 x 10~ 9 ). In that case, the positron signal from gravitino decay is indistinguishable from the 
background. 

have verified, is independent of the gaugino or gravitino masses used. In conclusion, the 
range of parameters that is compatible with gamma ray observations yields a negligible 
flux of positrons from gravitino decays. 

5.2 Antiprotons 

Gravitino decays into W*t and Z*v T produce antiprotons via gauge bosons fragmentation. 
We have used PYTHIA to obtain the predicted antiproton spectrum for several values of 
gaugino and gravitino masses. Figure [L4| shows the antiproton spectrum (before propaga- 
tion) for m^ = 70 GeV and different values of M\ . The number of antiprotons produced at 
a given energy is observed to increase with gaugino masses thanks to the larger branching 
ratio into the three-body final states. At 10 GeV, for instance, dNp/dE is about four times 
larger for M\ = 1 TeV than for M\ = 100 GeV. As before, this enhancement is saturated at 
M\ ~ 1 TeV, above which the antiproton flux no longer increases with M\. In contrast with 
the positron spectrum, the antiproton spectrum decreases at high energies rather quickly. 



The propagation of antiprotons in the Galaxy is described by a diffusion equation 
formally analogous to that for positrons, equation (|5.1|). Under certain conditions, see [21] 
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Figure 14: The antiproton spectrum (before propagation) as a function of the energy for different 
values of Mi . 



for details, the equation can be solved analytically to obtain 



d> p -(T,r) 



Vn 



An m 



v^ dm 



G 



dT 



(5.6) 



where i runs over the three-body final states W*t and Z*v T , T = E — nip is the antiproton 
kinetic energy, and Vp is the antiproton velocity. The function R(T) encodes all the as- 
trophysics and depends, in particular, on the propagation parameters. In our calculations, 
we will use the MED propagation model |ftq| , for which 5 = 0.70, Kq = 0.0112 kpc 2 /Myr, 
L = 4 kpc, and V conv = 12 km/s. It must be kept in mind, however, that the antiproton 
flux is very sensitive to the choice of propagation parameters. Indeed, an order of magni- 
tude variation with respect to the result for the MED model seems to be typical. For that 
reason, we will not try to derive constraints from antiproton data but simply illustrate how 
large the antiproton flux from gravitino decay can be. 

Before doing so, we must take into consideration that solar modulation effects can 



modify the observed antiproton flux at low energies. In the force field approximation [24] 
the antiproton flux at th< 
antiproton flux ($L ) by 



the antiproton flux at the top of the Earth's atmosphere ($J ) is related to the interstellar 



$ 



TOA 
P 



(T 



joA) = 2^ OA + TJ OA 
2m p Ti S + Tf OA 



(5.7) 



where Tjs, Ttoa are the antiproton kinetic energies respectively at the heliospheric bound- 
ary and at the top of the Earth's atmosphere. They are related as Tig = Ttoa + <Pf, with 
<j)p being the solar modulation parameter, which we take equal to 500 MV. 

In order to compare with available data, we must compute the antiproton to proton 
ratio rather than simply the antiproton flux. To do so, we need to know the proton 
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Figure 15: The antiproton to proton ratio as a function of the energy for different values of Mi. 
In this figure we set mp; — 70 GeV and the value of £ T is the maximum allowed by the gamma ray 
constraints (£ T = 7.9 x 10~ 10 ,2.4 x 10- 9 ,4.7 x 10~ 9 ,7.9 x 10~ 9 for Mi = 100,300,600, 1000 GeV, 
respectively). 



and antiproton backgrounds. The antiproton background, taken from p7| , 21 j, can be 
parametrized as 



<£bkg = 1 Q-1.64+0.071og 10 r-log? T-0.021og? T+0.0281ogf T 



whereas the proton background is given by [28] 

<f,bkg _ 6.88-0.828 log T-0.382 log 2 T+0.014 log 3 T+0.012 log 4 T-0.00097 log 5 T 



(5.8) 



(5.9) 



with T in GeV and the fluxes in m~ 2 s -1 str -1 GeV~ . The proton and antiproton fluxes 
are, With them and the predicted antiproton flux from gravitino decays, we can compute 
the expected antiproton-to-proton ratio at earth and compare it with the recent PAMELA 
measurements [29]. 



Figure 15 shows the p/p ratio as a function of the kinetic energy for mp = 70 GeV and 



different values of M\. For comparison, the PAMELA antiproton data and the predicted 
background are also shown. For each Mi the value of £ r used is the maximum that is 



allowed by the gamma ray constraint -see figure 1C. Notice that the different lines are in 



principle distinguishable from one another and from the background. For gravitino masses 



closer to Myy, the antiproton flux is even larger, as shown in figure 16. It is analogous 
to figure 15 but for rag = 80 GeV. Notice that, for Mi = 1 TeV, the antiproton-to- 
proton ratio deviates significantly from the PAMELA data points in both figures. For 
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Figure 16: The antiproton to proton ratio as a function of the energy for different values of M\. 
In this figure we set rriQ — 80 GeV and the value of £ T is the maximum allowed by the gamma ray 
constraints (£ T = 7.3 x 10~ 10 ,2.2 x 10~ 9 ,4.4 x 10~ 9 , 7.3 x 10~ 9 for M x = 100,300,600, 1000 GeV, 
respectively) . 



the MED propagation model we use, such models are clearly incompatible with present 
data. In the future, a better knowledge of the propagation parameters may allow us 
to use antiproton measurements to constrain the supersymmetric parameters directly. As 



indicated by figures |15| and |16|, the antiproton signal from gravitino decay can be significant 
even if the constraints from gamma rays are taken into account. 



6. Conclusions 

In the context of supersymmetric models with gravitino dark matter and bilinear R-parity 
violation, we have studied the indirect detection prospects of gravitinos for mg < My/. The 
main novelty in our analysis is the inclusion of gravitino decays into the three-body final 
states W*t and Z*v T , which were recently shown to modify in a significant way the gravitino 
lifetime and its branching ratios. As a result of these new contributions, the expected 
fluxes of gamma rays, positrons and antiprotons from gravitino decays are considerably 
altered. We first analyzed in detail the total gamma ray spectrum from gravitino decays, 
in particular the dependence of the new continuum contribution on gaugino and gravitino 
masses. Then, we computed the expected gamma ray flux from gravitino decays and derive 
new bounds on £ T and on the gravitino lifetime from recent FERMI data. The gravitino 
lifetime was shown to be constrained by such data to be larger than about 10 27 s to 10 29 s, 
depending on the gravitino and gaugino masses. We also investigated the indirect detection 
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of gravitinos via antimatter signals. In the range of gravitino masses we consider, they are 
entirely due to the three-body final states and had not been studied before. After analyzing 
the positron and antiproton spectra from gravitino decay and obtaining the expected fluxes 
at earth, we demonstrated that only the antiproton signal can be significant enough to be 
of interest for future experiments. 
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A. Analytical formulas for the three-body decays 

The gravitinos can decay to / + / and neutrino mediated by virtual photon or Z-boson, 
or / + /' and chargeed lepton by virtual VF-boson. Here we summarise the main formulas 
for the gravitino three body decays. 

The differential decay rate for the three body-decay process is given by 2 

dV N r 



dsdt 2567r 3 m ( =; 3 



|M| 2 , (A.l) 



where N c is the color factor (3 for q q final states and 1 for lepton pairs. The parameters 
s and t are the invariant masses of the // (or //') and fu systems respectively. We 
can define the masses and momenta of the relevant fields as {m^^p = k\ + k% + q) for 
Gravitino(Vv), (mi,ki) for a fermion (/), (m.2,&2) for the other fermion (/ or /'), (7713, q) 
for the neutrino or charged lepton and (my, k = k\ + k%) for an intermediate vector boson. 
Then the invariant masses s and t are 

s = (fci + k 2 ) 2 , t=(k 1 +q) 2 . (A.2) 

The range of these parameters are 

(mi + m,2) 2 < s < (mg — m 3 ) 2 , (A. 3) 



2 The analytic expression of the three-body decay are also shown in |p0|], |3ll] and 133 ] for R-parity 
conserving case. 
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and for fixed s, t are in the range 



tmin _; t _; T^maxt 



(A.4) 



where 



2 , _2 



tmax,min — ^3 1 Tlj T 



1 

27 



/ 2 2\/ 2 1 2\ 1 \l/2/ 2 2\\l/2/ 2 2 

(mg — s — m 3 ){s — m 2 + m^ ± A ' (s,mg , m 3 )A ' (s,m 1 ,m 2 

(A.5) 



The function A(x, y, z) is defined as 



A(x, y,z) = x + y +z — 2xy — 2yz — 2zx. 



(A.6) 



A.l G -)• 7*/Z* + v^ f + f + v 

The photon and/or Z-boson mediated three-body decay rate amplitude can be divided into 
three parts: photon- mediated, Z-boson mediated and the interference of the two, 



iM = iM~ t + iMz + iMi 



int- 



(A.7) 



For the range of the gravitino masses we consider, the photon-mediated diagram is subdom- 
inant and the amplitude is well approximated by the Z-mediated part. The two diagrams 
that contribute to gravitino decay into Z*u T (— >■ ^ , ffv T ) are shown in figure 17. 




c 




~("t) 



V t 



Figure 17: The two diagrams that contribute to gravitino decay into Z*v T 



The evaluation of those diagrams yields 

g 2 £ 



\Mi 



UMp 2 cos 2 9 W 3m 3 Mz 4 



G 



T z 2 M z 2 +(M z 2 -s) 2 



TrZ0 > (A.8) 



where g is SU(2) gauge coupling, 9\y is weak mixing angle and Y z is the decay width of 
Z-boson. The trace part is given by, in terms of the invariant masses s and t and neglecting 
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the masses of the final particles, 

TrZO 

= (30C A 2 sM z 6 + 30C v 2 sM z 6 - 36C A C v sM z 6 + 24C A 2 tM z 6 + 24C v 2 tM z G ) ,.. ( . 
+ (-30C A 2 s 2 M z G - 30CV VM Z 6 + 36C A C v s 2 M z 6 - 2AC A 2 t 2 M z 6 - 2AC v 2 t 2 M z t 






-48C A 2 stM z 6 - ±8C v 2 stM z 6 + 72C A C v stM z 6 ) 



m^, 



+ (2AC A 2 st 2 M z 6 + 24Cy 2 st 2 M z 6 + 24C A 2 s 2 tM z 6 + 24C7y VtM z 6 ) w,-, 



^ m^ 



+ E/-^ [(-38C A 2 sM z 5 - 38C v 2 sM z 5 + 76C A CysM z 5 , „ Y , 

+ (22C A 2 s 2 M z 5 + 22Cy 2 s 2 M z 5 - 92C A Cys 2 M z 5 + 14C A 2 siM z 5 + 14Cy 2 sfM z 5 



-l52C A C v stM z 5 ) m P A 



G 

5 



+ (WC A 2 s 3 M z 5 + 16Cy 2 s 3 M z 5 + 16C A Cys 3 M z 5 - 14C A 2 si 2 M z 5 - 14CV 2 si 2 M z 

i d 



-WC A 2 s 2 tM z 5 - WC v 2 s 2 tM z 5 + 32C A C v s 2 tM z 5 ) m^ 2 



-4C A 2 M z 5 s 2 t 2 - 4C v 2 M z 5 s 2 t 2 - AC A 2 M z 5 s 3 t - 4C v 2 M z 5 s 3 t] 
+ XJzz 2 [(UC A 2 sM z 4 + 14Cy 2 sM z 4 - 28C A C v sM z 4 ) m g 7 

+ (-10C A 2 s 2 M z 4 - 10CV 2 s 2 M z 4 + 20C A C v s 2 M z 4 - 2QC A 2 stM z A - 26C v 2 stM z 4 

+5QC A C v stM z 4 ) m d 5 
+ (2C A 2 s 3 M z 4 + 2CV 2 s 3 M z 4 + 20C A C v s 3 M z 4 + 26C A 2 st 2 M z 4 + 26C v 2 st 2 M z 4 



1 



+34C A ViM z 4 + 34CV VtM z 4 + 16C A C v s 2 tM z 4 ) m g 



+ (-6C A 2 s 4 M z 4 - 6Cy 2 s 4 M z 4 - 12C A Cys 4 M z 4 - 8C A Vi 2 M z 4 - 8Cy Vt 2 M z 4 
-8C A VtM z 4 - 8C v 2 s 3 tM z 4 - 24C A C v s 3 tM z 4 ] mg) . 

(A.9) 

Here U^ z is the zino-zino mixing parameter (see e.g. pi for its definition), which we 
assumed to be real. The coefficients Cy and C A are defined as 

Cv = \r 3 (f) - Q(f) sin 2 9w, C A = -~T 3 (f), (A.10) 

where T3 and Q are the third component of weak isospin and electric charge, respectively 
and we used Pl = (1 — 7s)/2. For example, for the left-handed neutrino, Cy = 1/4 and 
C A = -1/4. 

A.2 G ->• W* + T~ -)• / + f + T~ 

The two diagrams that contribute to the decay of the gravitino into W*r (Y^t ff' T ) are 
shown in figure [j]. The resulting squared amplitude for that process is given by 

W^V = 4i^i ^ ,r „ \ = —TrWO, (A . n) 



64A V 3m^M W A 



"G 



T 2 W M W 2 + (M w 2 - s) 
where g is SU(2) gauge coupling and Tyy is the decay width of VF-boson 
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The trace part is given by, neglecting the masses of final particles, 

TrWO 
= (96sM w 6 + 48tM w 6 )m d 5 + (-96s 2 M w 6 - 48t 2 M w 6 - 168stM w 6 ) mg 



+ (48st 2 M w 6 + 48s 2 tM w 6 )m d 

+ L7~~ [-l52M w 5 sm d 6 + (136s 2 M w 5 + l80stM w 5 )m d 4 

+(16s 3 M w 5 - 28st 2 M w 5 - 52s 2 tM w 5 ) m^ 2 - 8M w 5 s 2 t 2 - 8M w 5 s 3 t] 
+ £/~~ 2 [56M w *sm d 7 + (-40s 2 M w 4 - W8stM w 4 )m d 5 

+(-16s 3 M w A + 52st 2 M w 4 + 52s 2 tM w A ) m g 3 + (8M w 4 sH - 16M w 4 s 2 t 2 )m 5 ] , 

(A.12) 

where U^^, the wino-wino mixing parameter (see e.g. J?]]), was assumed to be real. 
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